Properties of edge states in spin-triplet two-band superconductor 
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Motivated by Sr2Ru04 the magnetic properties of edge states in a two-band spin-triplet supercon- 
ductor with electron- and hole-like Fermi surfaces are investigated assuming chiral p-wave pairing 
symmetry. The two bands correspond to the a-/3-bands of Sr2Ru04 and are modeled within a tight- 
binding model including inter-orbital hybridization and spin-orbit coupling effects. Including super- 
conductivity the quasiparticle spectrum is determined by means of a self-consistent Bogolyubov-de 
Gennes calculation. While a full quasiparticle excitation gap appears in the bulk, gapless states 
form at the edges which produce spontaneous spin and/or charge currents. The spin current is 
the result of the specific band structure while the charge current originates from the superconduct- 
ing condensate. Together they induce a small spin polarization at the edge. Furthermore onsite 
Coulomb repulsion is included to show that the edge states are unstable against the formation of a 
Stoner-like spin polarization of the edge states. Through spin-orbit coupling the current- and the 
correlation-induced magnetism are coupled to the orientation of the chirality of the superconducting 
condensate. We speculate that this type of phenomenon could yield a compensation of the magnetic 
fields induced by currents and also explain the negative result in the recent experimental search for 
chiral edge currents. 
74.70.Pq,: 



I. INTRODUCTION 

Among the unconventional superconductors Sr2Ru04 
has received special interest for its unique and most in- 
triguing properties which resemble in some aspects the 
spin-triplet superfluid '^He [H, Despite the relatively 
low transition temperature of Tc ~ 1.5 K and the fragility 
of the superconducting phase against disorder, a large 
bulk of experimental data is available nowadays which 
gives strong evidence that a pairing state with so-called 
chiral p-wave symmetry is realized [1, 0]. This state 
breaks time reversal symmetry implying magnetic prop- 
erties due an intrinsic angular momentum of the Cooper 
pairs. Internal magnetism of the superconducting phase 
has indeed been observed by zero-field relaxation of muon 
spins As an odd-parity spin-triplet pairing state the 
gap function of the chiral p-wave state can be represented 
in the vector representation 



d{k) = Aoz(/Ca; ± '^ky), 



(1) 



where the z-axis orientation of d indicates an equal-spin 
pairing state within the x-y-planc. This is the analog of 
the A-phase of "^He and is two-fold degenerate, i.e. it has 
positive or negative chirality depending on the sign of the 



* I imai@phy.saitama-u.ac.jp] 



angular momentum L^, where = ±1 for the orbital 
dependences kx ± ifcy of the gap function. 

The chirality of this pairing state is a topological prop- 
erty and generally yields edge states with chiral proper- 
ties ^ . These states are so-called Andreev bound states 
localized near the surface of a sample penetrating on the 
length scale of the coherence length |7]] . These currents 
carry a finite supercurrent whose direction in the x-y- 
plane is connected with the sign of L^. It has long been 
suggested that such orbital current might be detectable 
by sensitive local probes [1, H, 0| • So far all attempts to 
find the magnetic fields induced by these currents, us- 
ing scanning Hall probes and SQUID microscopes, have 
led to negative results 10, 11]. On the other hand, the 
presence of surface subgap states has been demonstrated 
by quasiparticle tunneling spectroscopy [l^.[l3j. Thus, it 
remains an unresolved puzzle that the predicted currents 
lie below the detectable limit, in contradiction with the 
theoretical estimates [l^ . 

It has been early realized that superconductivity in 
Sr2Ru04 is more complex due to multi-orbital band 
structure 15]. Indeed it has been suggested that due 
to the different character of the bands incorporated in 
this material also the properties of the surface states 
could be influenced [11]. The three bands crossing the 
Fermi energy derive from the the 4d-i2g-orbitals of the 
Ru ions. These orbitals {dyz,dzx,dxy), disperse via tt- 
hybridization with intermediate O 2p-orbitals giving rise 
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to the characteristic two-dimensional band structures 
which can be easily derived from the corresponding tight- 
binding model in the x-y-plane of the layered crystal 
structure, invoking very weak dispersion along the z-axis. 
The three resulting bands separate into two groups the 
Q!-/3-bands belonging to the dyz- and dzx-oihitals with 
negative parity under reflection at a;-y-plane, in contrast 
to the dxy-oihital which yields the 7-band. Several exper- 
imental results indicate that the 7-band with its electron- 
like Fermi surface is dominant for superconductivity 
while the other two bands are responsible for the strong 
incommensurate magnetic correlation due to their nearly 
one-dimensional character [13, HI]- Raghu et al. have 
pointed out that in terms of topology the hole-like a- 
and the electron-like /3-band compensate each other such 
that the Chern numbers of the topology introduced by 
the chiral p-wave state adds up to zero (16!]. Based on 
this finding they give an argument for why the surface 
currents would be unobservable, if these two bands were 
the dominantly superconducting ones. 

In the present study we concentrate also on the a-/3- 
bands. Our aim is to show how superconducting and 
magnetic properties of these two bands could combine in 
a peculiar way to give rise to novel edge states in the 
superconducting phase. It was noticed earlier that even 
the normal state has unusual properties at surfaces with 
a normal vector in the x-y-plane. Spin-orbit coupling 
combined with the inplane hybridization structure of the 
dyz-dzx-OThitals gives rise to surface spin currents p^ . 
which are intimately connected with the presence of a 
strong anomalous Hall effect proposed for Sr2Ru04 (20I - 
[2^ . Interestingly, together with the chiral edge current 
this spin current can generate a finite spin magnetiza- 
tion at the edge, which is related to the spin Hall ef- 
fect. This magnetization may be enhanced due to strong 
magnetic correlations intrinsic to the a-/3-bands. In our 
study we discuss a simple two-band model featuring the 
Q!-/3-bands and incorporating both superconductivity and 
magnetism. The 7-band may support and induce su- 
perconductivity on the Q!-/3-bands but is ignored in the 
present treatment for simplicity. The spin magnetism oc- 
curring at the edge contributes to the magnetization at 
the surface and may, in principle, act as a compensation 
to the orbital supercurrents to diminish the signal to be 
detected by scanning magnetic probes. 

This paper is organized as follows. We construct the 
effective model in the next section. In Sec. IIIIl we show 
the spin and charge current. In addition to the currents, 
the band structure and magnetization in the supercon- 
ducting phase are shown in Sec. |lVl Summary and dis- 
cussions are given in Sec. |Vl 



II. MODEL 

In this section we introduce the basic model which we 
investigate for its edge properties. Edges can be most 
easily addressed in a ribbon-shaped system, an infinite 
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FIG. 1. Lattice structure of i-leg ribbon with two px and Py 
orbitals. t {t') stands for the hopping integral between same 
(different) orbitals. 



ladder with L legs as shown in Fig|TJ providing two edges. 
The two orbitals can be viewed like Px- and Pj,-orbitals 
(FigH]) representing dzx- and dj^x-orbitals, respectively. 
Obviously, the nearest neighbor hopping is only possi- 
ble via intra-orbital hybridization with a hopping matrix 
element t. This yields for both orbitals one-dimensional 
dispersion along the x- and j/-direction, respectively. The 
inter-orbital hybridization involved next-nearest neigh- 
bor hoping of strength (hopping matrix element t') as 
depicted in FigHJ 

Including onsite spin-orbit coupling and interactions 
the corresponding Hamiltonian is written as 



H — -ffdd + Hdd' + Hso + -ff a + Hi- 

with 



(2) 



L-l ' \ 

+ 4lpy,y(=WPya + h.C. 1 - M XI 

l—l / ilma 

i,m.(7 \ l — l 

- Yl 4m^C,+i/_irna + h.C. 1 , 
1^2 ) 



ia I 

Hr = Ur ^ num^nurni + ^ ^^z/mt^i/m^ 
ihn ihn 
-\-{K-p — t/j-) ^ ^ ^ilma'^ilrfKT 1 
ilm{<fh)(T 

(3) 

where cj^^^ {cumcr) is a creation (annihilation) operator 
for electron as site (i, /), orbital m = {px,Py) and spin a 
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where i represents the coordinate along the a:-direction 
and I the leg index and m = Py,Px)- Moreover, /i and 
A denote the chemical potential and the spin-orbit cou- 
pling strength. Several interaction terms are introduced. 
The effective attractive nearest neighbor interaction term 
with Ua (< 0) will be used to form p-wave Cooper pairs. 
The repulsive onsite interaction among the electrons in- 
cludes intra-orbital Ur, inter-orbital Kr and Hund's rule 
coupling Jr. For the latter we ignored the exchange and 
pair hopping part in order to simplify the discussion with- 
out changing the qualitative outcome. We furthermore 
assume the standard relation Ur — Kr + 2 Jr. These 
repulsive terms are important for the magnetic correla- 
tions. 

Neglecting the interaction terms we can determine the 
band structure leading to the Fermi surfaces shown in 
Fig 12] where the dashed lines indicate the purely one- 
dimensional Fermi surfaces in the absence of inter-orbital 
hybridization and spin-orbit coupling (t' = A = 0). We 
have adjusted here the chemical potential as to obtain 
a particle density of 8/3, close to the experimentally es- 
timated band filling. Obviously even with intcr-orbital 
coupling the nesting property of the Fermi surfaces pre- 
vails, while the electron-like /3 Fermi surface centered at 
the F-point and the hole-like a Fermi surface centered at 
the M-point are formed. 

In the following the interaction terms will be decou- 
pled by mean fields. For the attractive nearest-neighbor 
interaction we use the BCS-type mean field with the su- 
perconducting order parameters defined as 



ma.m' a' 



\Ci-\-ll'm(7Cil'm' a' 



(4) 
(5) 



Note these mean fields are located on the bonds 
rather than on the sites. We will restrict here 
on the spin-triplet channel with inplane equal-spin 
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FIG. 2. Fermi surfaces of two-dimensional bulk system for 
t' = O.lt and A = O.lt (solid line) and t' = A = (dashed line), 
in which particle number n = 8/3 and Ua = Ur ~ Jr = 0. 



pairing using only the intra-orbital mean fields 

^Lma,m<j ^ Srn,p^{Ci+llp,aCilp^a) = and Af^^^^^- = 

Sm,Py{cii+ipyaCiipya) = Af whilc ah other pairing decou- 
plings are assumed to vanish. This leads to the chiral 
p-wave phase as the most stable superconducting phase. 
Since the definition of is not symmetric with respect 
to L/2, the redefinition of the gap function is given by 



Af = 



Af/2 



(At 



Ar/2 



{I 
{I 



1) 
L) 



Af)/2 (otherwise) 



(6) 



For the onsite interactions we apply the usual Hartree- 
Fock approximation with particle density and the spin 
polarization as follows. 



where 



nir 



{nur 



(4 



(7) 
(8) 



(9) 



Our ribbon geometry allows us to consider the system 
as homogeneous along the x-direction, but implies the 
spatial dependence of these mean fields in the transverse 
direction. 



III. SPIN CURRENTS IN THE NORMAL STATE 

Before discussing the superconducting phase we first 
address the surface properties of the normal state for 
Ua = 0. To illustrate the behavior of the system qualita- 
tively we use mainly the following model parameters: the 
ribbon L = 100, the inter-orbital hopping matrix element 
t' — O.lt, and the electron density n = 8/3. 

The topology of intra- and inter-orbital hybridization 
and spin-orbit coupling leads to spin currents at the sur- 
face. This can be observed on the level of the single- 
particle part of the Hamiltonian so that we neglect the 
interaction terms. We introduce the momentum k along 
the ribbon and define the new electron operators 




-ikxi 



(10) 



with Xi being the x-coordinate of site (i,Z) (i^; the 
length of the ribbon assuming periodic boundary con- 
ditions along x-direction). 

The spin dependent current operator along the a;-axis 
consists of two parts. 



_■(!), -(2) 



(11) 



corresponding to an intra-orbital contribution through 
nearest-neighbor hopping along the a;-axis. 



(1) 



klp^a'^klpaicr, 



(12) 
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FIG. 3. Spin-dependent current in normal phase for L — 100, 
t' = \ ^ O.lt, [/r = Jr = and n = 8/3. Solid (dashed) line 
stands for spin up (down) current. The inset shows magnifi- 
cation near the edge. 

and an inter-orbital part with next-nearest neighbor hop- 
ping diagonal in the square plaquettes, 

k 

). (13) 

The expression given here for is not symmetric 
with respect to I — L/2. Therefore we redefine it as 

, f jL'V2 a = l) 

il'^ = if 1./2 (l^L) . (14) 

[ (jf\.+jf')/2 (otherwise) 
so that the spin-dependent current operator is now given 

by 

Ji.-Jl^+Jlf, (15) 

which does not affect the results qualitatively. The 
charge and spin currents are then given by 

^- = 4E{3i^^+ir}' (16) 

It is straightforward to determine the groundstate of 
the normal phase and to calculate the expectation values 
for the currents. The spin-dependent current is depicted 
in Fig [H As time reversal symmetry is conserved in the 
normal phase there is no charge current running. On the 
other hand, we find a spin current at the edges which 
extends towards to center of the ribbon with a charac- 
teristic oscillation. 

The presence of such spin currents can be anticipated 
by rather simple considerations in the topology of the 
single-particle Hamiltonian as is revealed, if we exam- 
ine the phase structure of a square plaquette, as drawn 
in FiglH Separating a single triangle within the plaque- 
tte, we see that in the combination of nearest-neighbor 



intra- and next-nearest neighbor inter-orbital hopping 
with the onsite spin-orbit coupling a single particle of 
spin a picks up a spin-dependent phase (j)^ = cr7r/2 mov- 
ing once around the triangle(cr being the spin ±1 with 
quantization axis z). It is straightforward to see that 
each single particle state carries, in general, a finite spin- 
dependent currents of equal magnitude around each tri- 
angle, but opposite sign for the two spins cr. This gives 
rise to a finite spin but a vanishing charge current. As 
they have the same orientation for all four triangles in- 
scribed into the square, the diagonal bonds have no net 
current, while the bonds on the square carry a spin cur- 
rent all around. Taking the whole system, these spin 
currents cancel out in the bulk, but remain at the edges 
due to lack of compensation. 

Alternatively, the phase winding on the triangles can 
be represented as a "spin-dependent flux" density 
with —$4, due to time reversal symmetry. We can 

define a spin dependent current by 

j,{r) = V X z$,(r) (18) 

which is yields by symmetry j'^ = j-^ -\- ji = and, in 
general, a non- vanishing spin current density, if $(j(t") is 
not uniform. The latter condition is satisfied at the edges 
of our system. 

The oscillations in magnitude of the spin current par- 
allel to the edge can be interpreted as Friedel oscillations 
as they correspond roughly to the nesting vector in the 
band structure, which also is the wave vector of the dom- 
inant spin correlation. 



IV. SUPERCONDUCTING PHASE 

In the following we treat the system within the super- 
conducting phase based on the model given by Eas. (|2l3l) . 
For this purpose we will solve the corresponding mean- 
field in the ribbon with spatial resolution self-consistently 
for Ua — — 1.5t at zero-temperature. The other param- 
eters are as in Sec. IIIII We choose a strong attractive 
potential Ua in order to obtain a rather small coherence 
length of a few lattice constants for computational pur- 
pose to keep the system size small. 
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FIG. 4. Spin-dependent chirality. 
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FIG. 5. (Color online) Order parameters of superconducting 
state as a function of leg index / for X = O.lt and Ur ~ Jr ~ 0. 
The black (red) lines stand for Ur = 0.0 {Ur = 0.5t). The 
inset shows magnification near the edge. 



A. Superconducting order parameter and energy 
dispersion 

First we analyze the behavior of the superconduct- 
ing order parameter assuming spin-triplet pairing with 
inplane equal-spin pairing. In this case the most sta- 
ble pairing state has chiral symmetry in order to avoid 
nodes in the quasiparticle gap. The order parameter 
components and have a relative phase tt/2 (or 
the degenerate — 7r/2) and are of equal magnitude in the 
bulk, e.g. ReA^=ImA^ , as expected for the chiral p- 
wave state d — z{kx + iky). In the ribbon the order 
parameter components show a characteristic spatial de- 
pendence as depicted in FiglS] where ReA^ and ImA^ 
are shown as a function of the leg index I. The behav- 
ior at the edges originates from surface scattering effects, 
where A^ (even under reflection at the edge) slightly in- 
creases, while A^ (odd under reflection at the edge) is 
suppressed [1, Q . The pairing state relies on the attrac- 
tive nearest-neighbor interaction iJai while the repulsive 
onsite interactions in affects the order parameter very 
weakly. Intentionally we keep at this level these onsite 
interactions weak as they trigger a magnetic instability 
which we will discuss below. At this point a few com- 
ments are in order. Within our treatment the pairing 
state with ReA^=— ImA*' is degenerate in energy and 
corresponds to the time reversed state d — z{kx — iky). 
Spin-orbit interaction does not affect this general fact, 
as time reversal symmetry is broken spontaneously here. 
Hereafter we discuss mainly the d = z{kx + iky) state 
unless otherwise noted. 

Figure [5] shows the energy spectrum for Ur — Jr — ^ 
and Ua = — 1.25i and — 1.5t. The result shows a gapped 
"bulk" spectrum as expected for the chiral p- wave state. 
In addition two branches of subgap gap states appear 
corresponding to Andreev bound states of the two ribbon 
edges. Especially in Fig. [S] (b) it can be seen that these 
edge states are decoupled from the bulk spectrum. As the 




FIG. 6. (Color online) Energy dispersions near low-energy 
region; (a) Ua = — 1.25i and (b) Ua = — 1.5f for Ur = Jr ~ 0. 
Red (blue) line stands for edge state from I — 1 (I = L). 



two Fermi surfaces have electron- and hole-like character 
their topological Chern numbers cancel [l^ US] • There- 
fore each of the two edges contributes two zero-energy 
crossings of opposite chirality, such that these edge states 
are not topologically protected. 



B. Spin and charge currents 

We now analyze the equilibrium currents in the ab- 
sence of onsite repulsion. The spin dependent current 
densities along x-directions are shown in Fig 17] as a func- 
tion of the leg index I for both chiral p-wave states 
d = z(fca; ± ifcy) as . Note that since the current 
from off-diagonal (superconducting) term becomes less 
than one-tenth of jia, the contribution is negligible. 

Obviously the currents are associated with the edges. 
It is important to see that the main contribution to the 
currents originates from the intra-orbital (j)x) component 
jj'j;'* which exceeds clearly the inter-orbital as shown 
in the inset of FigIT] (for d ~ z{kx + iky)). 

Various symmetry properties can be stated for the cur- 
rents. The relation between positive and negative chiral- 
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FIG. 7. (Color online) Spin-dependent currents as a function 
of / for A = O.lt and Ur = Jr = 0. The black (red) lines 
stand for d = z{kx + iky) {d — z(kx — iky)) state and the 
solid (dashed) lines stand for spin up (down) contributions. 
The inset shows each component of the spin up current for 
d = z{kx + iky) state. 



ity is reflected by 



,■(+) _ _A-) 

Jlcr ~ -lis 



(19) 



The sign of the spin-orbit coupling constant A yields the 
relation, 

(signA = +1) - (signA = -1). (20) 



We find finite charge and spin current densities which 
is shown in FiglSl Obviously, the sign of the charge (spin) 
current depends on the sign of chirality of the supercon- 
ducting phase (spin-orbit coupling). 

Figure [9] shows the sums of spin and charge currents 
up to / = L/2, which is defined as 



L/2 

js(c) _ f{c) 
"^sum ~ / , "^l 

1=1 



(21) 



FIG. 9. Sum of spin current for various A. The inset shows 
the sum of charge current for Ur = Jr ~ 0. 



When the number of legs is sufficiently large, Jsum is 
saturated with increasing I from the edges, which corre- 
sponds to the total current at each edge. In comparison 
with the result of the charge current, the sum of spin cur- 
rent strongly depends on the amplitude of the spin-orbit 
interaction. 

Next, the effect of repulsive interactions on the spin 
and charge currents is discussed. Figure [10] shows the 
spin and charge currents with/without repulsive interac- 
tions. Although the amplitude of each current is slightly 
modified by repulsive interactions, the sign of the current 
at each / never changes. In particular, the spin current is 
almost independent of the amplitude of repulsive inter- 
actions, in which the contribution of the horizontal com- 
ponent (if^') is also dominant in the various repulsive 
interaction region. We stress that the currents near the 
edges fiow mainly between orbitals along the edges and 
the contribution of repulsive interactions to Px orbitals is 
rather inconspicuous. 



0.02 I . , . , . , . , r 




FIG. 8. Spin and charge currents as a function of I for A = O.lt 

and Ur = Jr = 0. The inset shows the flow direction of the FIG. 10. (Color online) Spin and charge currents for 
chiral charge current near the edges. with/without repulsive interactions near an edge for A = O.lt. 
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FIG. 11. (Color online) Spin polarization for various choices 
of A as a function of I. (a) and (b) stand for mip^ and mip^, 
respectively. Ur = 0-5t and Jr = 0. The insets show the 
magnifications near edge. 



C. Magnetic Properties at edges 



both spin-orbit coupling and the next-nearest neighbor 
hopping t' . 

This magnetic instability for the pj^-orbital is of Stoner- 
like nature, introducing a spin splitting of the surface 
states due to the exchange interaction. This aspect is 
clearly visible in Fig |12l which shows the spin polarization 
for different values of Ur, showing a stronger magnetiza- 
tion for larger interaction Ur- In sharp contrast, mip^ 
is almost independent of Ur- Consequently we may con- 
clude that the magnetizations carried by the two orbitals 
are of different origin. Through the system geometry 
and hopping topology the current-induced spin polariza- 
tion is dominantly associated with the Px orbital, while 
the Stoner-induced spin polarization can, obviously, be 
attributed to the py-orbitals. 

In the absence of spin-orbit coupling the magnetic or- 
der of the py-orbitals at the edges has a continuous de- 
generacy due to SU (2)-spm rotation. Spin-orbit coupling 
pins the orientation of the spin polarization to z-axis. 
This is connected with the fact that spin correlations 
involving the nesting wave vector yields a susceptibility 
favoring the z-axis polarization [Tsj . Thus, the remaining 
two-fold degeneracy (magnetization parallel and antipar- 
allel to z-axis) can now be broken by the chirality of 
the superconducting phase. The orientation of the spin 
polarization mip^ of p^j-orbital is determined by the di- 
rection of the chiral edge current, because mip^ is the 
result of the combination of spin and chiral edge current. 
Through spin-orbit coupling, inter-orbital hybridization 
and Hund's coupling between the two orbital, the spin 
polarization m/p^ generates a bias for the orientation of 
mipy . Thus, the spin magnetization of both orbitals is 
tied to the chirality of the superconducting phase. 



The essentially flat electronic bands of the edge states 
which cover a considerable range of momenta k along 
the ribbon (see FiglHl) derive mainly from the py-orbitals. 
This special role of the pj^-orbital is also compatible with 
the suppression of the -component of the supercon- 
ducting order parameter near the edges which goes hand 
in hand with the formation of subgap edge states. These 
edge states contribute a large density of states at zero 
energy at the edges and are, therefore, very susceptible 
to a "Fermi surface" instability. Indeed, turning on the 
repulsive interaction Ur yields a finite spin magnetization 
at the edges. The magnetization is uniform along the rib- 
bon and extends towards the interior with an oscillation 
characteristic to the nesting vector of the band structure 
(see FigHUandHll). 

Figure [TT] shows the spin polarization for various 
choices of the spin-orbit interaction for finite repulsive 
interaction Ur- The spin polarization oi orbitals for 
A = is very small, but increases with growing spin- 
orbit coupling A near the edges. On the other hand, the 
spin-orbit interaction does not affect the spin polarization 
mip much. Note, that the two moments uiip^ and mip 
associated with the two orbitals are connected through 



D. Induced magnetic fields 

Both the edge charge supercurrent and the spin mag- 
netization generate a net magnetic field. Let us here 
compare their magnitudes. Using the Maxwell's equa- 
tion V X B = fioj we find that magnetic field of the 
current is given by 

Btil) = MoE('^^') = -^^oE(Jp), (22) 
I' I' 

where a is the inplane lattice constant, is the mag- 
netic constant and (J^) is the dimensionless current den- 
sity along the x-direction on layer /. The result is shown 
in FigHni Note, that the screening effects have been ig- 
nored here such that the magnetic field is constant in- 
side the ribbon. We are only interested in a rough es- 
timate of magnitude of the field and Meissner screening 
would introduce counter currents suppressing this field 
on a length scale of London penetration depth. 

The spin polarization produces a local magnetization 
to the magnetic field, 

5I(0 = -^A^bK-"4), (23) 
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FIG. 12. (Color online) Spin polarization for several choices 
of Ur- (a) and (b) stand for niip^ and mip^, respectively. 
A = O.lt and Jr = 0. Both insets show the magnifications 
near the edge. 



where fiB is Bohr magneton and c the lattice constant in 
z-direction. For the mutual comparison of the two fields 
it is useful to consider the dimensional prefactors. Thus, 
we consider 



ha I 



a?c m* 



0(1) 



(24) 



taking fJ-B = he /2m and m* as the effective mass. 

In FigllS] the two magnetic fields can be compared 
both given in the corresponding units, fiQct/ha for 
and fio^B/a^c for B^. We observed that both contribu- 
tions are of similar magnitude and have opposite sign. 
Since the two units are of similar order the two fields 
tend to cancel each other. Note that this compensation 
occurs for both types of chiral domains, since the ori- 
entation of the spin magnetization is through spin-orbit 
coupling coupled with the chirality. 



V. SUMMARY AND DISCUSSIONS 

Among the three important bands of Sr2Ru04 two, the 
a- and /3-bands derived from the dyz and dzjj-orbitals, 
show remarkable magnetic properties. The special struc- 
ture of spin-orbit coupling induces here features of an 




FIG. 13. (Color online) Sum of spontaneous magnetic field for 
several choices of repulsive interaction. Solid (dashed) lines 
stand for the magnetic field from the charge current (spin po- 
larization) (A — O.lt and J,- = 0). Note no screening currents 
have been included in this calculation. 



anomalous Hall effect. In particular, the edges carry 
spin currents. We have analyzed possible consequences 
of these features on the superconducting phase. Assum- 
ing chiral p-wave symmetry for the Cooper pairs Andreev 
bound states at the surface not only induce surface super- 
currents, but also introduce a spin polarization. There 
are two components to the spin magnetism here. One is 
induced by spin Hall effect like features which give rise to 
a spin density imbalance at edges when simultaneously 
charge and spin currents are present. The other compo- 
nent nucleates due to the strong magnetic correlation in 
the a- and /3-bands of incommensurate spin density wave 
type, which benefit from the modification of quasiparticle 
spectrum at the surfaces. 

The combination of both shows that the chirality of the 
superconducting phase and the orientation of the spin po- 
larizations are coupled. The sign of spin-orbit coupling 
is so as to lead to opposite magnetizations for the su- 
percurrents and spin magnetic moments. We speculate 
here that this phenomenon might provide a way to reduce 
the surface magnetization of the superconducting phase 
through a compensation effect. In this way the negative 
result of the attempts to detect the magnetic fields of a 
chiral p-wave phase by means of local magnetic probes 
could be explained. One could hope that the observation 
of the surface spin polarization would be possible by a 
spin resolved tunneling experiment, which would provide 
a direct test for our scenario. 

Our model has certainly short-comings as the 7-band 
dominating most likely superconductivity in Sr2Ru04 
has been neglected here. Certainly the coupling between 
the two subsets of bands transfers superconductivity from 
the 7- to the a-/3-bands and, through spin-orbit and 
Hund's rule coupling, the spin-polarization is induced in 
the 7-band. In any case we expect a reduction of the bare 
surface current induced local magnetic fields. A further 
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issue is the effect of disorder on the discussed effects. 
Chiral edge states of the a-/3-hands are not topologi cally 
protected, since the total Chern number vanishes [16| . 
Thus, the situation may be ahered by surface roughness 
and order disorder effects. These questions we will be ad- 
dressed in future studies. We would like to mention also 
that recently, the aspect of the anomalous Hall effect in- 
trinsic to the a-/3-band has recently also been discussed 
by Taylor and Kallin (2^ in the context of the polar Kerr 
effect m. 
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